Domain growth morphology in curvature driven two dimensional coarsening 
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We study the distribution of domain areas, areas enclosed by domain boundaries ("hulls"), and 
perimeters for curvature-driven two-dimensional coarsening, employing a combination of exact anal- 
ysis and numerical studies, for various initial conditions. We show that the number of hulls per 
unit area, nh{A, t) dA, with enclosed area in the interval [A, A -f dA), is described, for a disordered 
initial condition, by the scaling function nh{A,t) — 2ch/(A + Xht)^, where Ch ~ l/87r\/3 ~ 0.023 
is a universal constant and is a material parameter. For a critical initial condition, the same 
form is obtained, with the same Ah but with replaced by C(,/2. For the distribution of domain 
areas, we argue that the corresponding scaling function has, for random initial conditions, the form 
nd{A,t) = 2cd{Xdty ~^ /{A + Xdty , where Cd and Ad are numerically very close to Ch and Ah re- 
spectively, and r' — 187/91 ~ 2.055. For critical initial conditions, one replaces Cd by Cti/2 and the 
exponent is r = 379/187 « 2.027. These results are extended to describe the number density of 
the length of hulls and domain walls surrounding connected clusters of aligned spins. These predic- 
tions are supported by extensive numerical simulations. We also study numerically the geometric 
properties of the boundaries and areas. 



I. INTRODUCTION 

In this work we obtain some exact and new results 
for the coarsening dynamics [l| of a non-conserved scalar 
field in two dimensions, demonstrating, en passant, the 
validity of the dynamical scaling hypothesis for this sys- 
tem. We study the morphology of the domain structure, 
which is illustrated in Fig. [T] for the coarsening of the 
two-dimensional Ising model (2c?IM) on a square lattice 
quenched from an equilibrium state at Tq > Tc- 

A domain is a region of connected aligned spins. Each 
domain has one external perimeter which is called the 
hull. The hull enclosed area is the total area within this 
perimeter, i.e. the domain area plus the area of any in- 
ternal sub-domain. The domain perimeter is the total 
length of the interface between the chosen domain and 
the neighbouring ones - including the hull and internal 
borders. See Fig. [2] for a sketch explaining these defini- 
tions. 

The paper deals primarily with the distributions of 
two characteristic areas, the domain area and the hull 
enclosed area, and their associated lengths, the domain 
wall perimeter and the hull length. 

Naively, one may imagine that coarsening is basically 
due to the coalescence of small domains that form larger 
ones. However, in two dimensional curvature driven 
coarsening coalescence processes are quite unimportant 
as shown by the AUen-Cahn result. All the interfaces 
move with a local velocity that is proportional to the 
local curvature and points in the direction of decreas- 
ing the curvature; therefore, interfaces tend to disappear 
independently of one another. This is the reason why 
we first focused on the statistics of hull enclosed areas, 
quantities that depend on the motion of a single and con- 
nected interface, and not on the statistics of the more 
natural domain areas. Next we expressed the statistics 
of the domain areas in terms of the simpler and more 



clear statistics of hull enclosed areas. 

Hull enclosed and domain areas have distributions 
that, at late times after the quench into the ordered 
phase, exhibit, according to the scaling hypothesis, the 
scaling form n{A, t) = t^"^ f{A/t), where n{A, t)dA is the 
number of hulls (domains) per unit area with area in the 
range [A, A + dA). The argument of the scaling function 
arises from the fact that the characteristic length scale 
is known to grow as t^/"^, so the characteristic area (of 
hulls and domains) grows as t. The scaling function f{x) 
will be different for domains and hulls. The prefactor t~'^ 
follows from the fact that there is of order one hull (or do- 
main) per scale area. Our analytical result is an elegant 
application of the Gauss-Bonnet theorem and, it should 
be emphasized, its simplicity relies on the dimensionality 
of the system being 2. Indeed, in three dimensions, the 
time variation of the hull enclosed volume depends on a 
characteristic size of the domain. A similar dependence 
on the dimension is also observed in the von Neumann's 
law for cellular systems, whose simple form, independent 
of any linear size of the system, is also only observed in 
two dimensions 0, 

In this paper we derive these scaling forms from first 
principles (i.e. without recourse to the scaling hypothe- 
sis), and determine explicitly the scaling functions. Some 
of our results have appeared earlier in Letter form . 

Hulls and domain boundary lengths are themselves dis- 
tributed quantities related in a non-trivial manner to 
their corresponding areas. In this paper we examine the 
geometry of these structures and we derive the number 
density of hull and domain wall lengths showing that 
these distributions also satisfy scaling. 

The organisation of the paper is as follows. In Sect. |TT] 
we recall known results about the equilibrium distribu- 
tion of hull enclosed and domain areas at critical percola- 
tion and critical Ising initial conditions. We also summa- 
rize known results about the equilibrium distribution of 
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domain walls and hulls and their geometrical relation to 
their associated areas. In Sect. lIIII we derive some generic 
results that stem from a number of sum rules and the 
use of the scaling hypothesis. In Sect. II VI we explain the 
analytic derivation of the time-dependent hull enclosed 
and domain area distributions. These arguments do not 
rely on any scaling hypothesis but rather demonstrate 
its validity. In Sect. |V] we show our numerical results for 
the statistics of areas in the 2dlM evolving with Monte 
Carlo dynamics. Section rvTI is devoted to the analysis, 
both analytical and numerical, of the geometry of hulls 
and domain walls during the dynamics and their relation 
to their corresponding areas. Finally, in the Conclusion 
we discuss future studies along these lines. We also add 
two appendices in which we describe the algorithm used 
to identify and count the hull enclosed areas, and for the 
sake of comparison we present the distribution of domain 
lengths in the onc-dimcnsional Ising model. 




The equilibrium hull enclosed area distributions at perco- 
lating criticality and Ising criticality, have been computed 
by Cardy and Ziff in two dimensions: 

(A n^ / 2c/i/A^ , critical percolation, , . 
"^^^'"^^ \c^M2 , critical Ising. 

These results are valid for ^ A ^ i^, with Ao 
a microscopic area and the system size. Note also 
that we are taking an extra factor 2 arising from the 
fact that there are two types of hull enclosed areas, cor- 
responding to the two phases, while the Cardy-Ziff re- 
sults accounts only for clusters of occupied sites (and not 
clusters of unoccupied sites). nh{A^O)dA is the num- 
ber density of hulls per unit area with enclosed area in 
the interval {A, A + dA) (we keep the notation to be 
used later and set t = 0). The adimensional constant 
Ch is a universal quantity that takes a very small value: 
Ch = 1/(877^/3) « 0.022972. The smallness of c,, plays an 
important role in the analysis of Sect. IV. 




FIG. 1: Snapshots of the 2d Ising model at time t = 32 MCs 
after a quench from infinite temperature. To 00, to the 
working temperature T = 1.5 ~ 0.66Tc. In the left panel 
we show the raw data, where the domain structure as well as 
the thermal fluctuations within the domains are visible. In 
the right panel we show the boundaries between regions of 
opposite sign in the configuration to the left, using a variant 
of the algorithm in Ref. [sj to make the domain structure 
clearer. 



FIG. 2: (Colour online.) A sketch of a configuration with 
two concentric and circular interfaces with radius _Ri and R2 
is shown to illustrate the definition of hull enclosed and do- 
main areas as well as hull and domain-wall perimeters. This 
configuration has two hull enclosed areas A'j^'> = ttRI and 
A^f^^ = ttR2 and two domains with areas A^J''' = ■k{Ri — R2) 
and A'^ ' = TrRi. There are two hulls with length pf^ — 2-kR\ 
and p^^' = 2nR2, and two domain walls with length = 
27r(i?i -I- R2) and p^^' = 2-kR2. 



II. EQUILIBRIUM DISTRIBUTIONS 

In this Section we summarize known results on the 
equilibrium hull enclosed and domain area distributions 
as well as the number density of perimeters, both at 
critical percolation and critical Ising conditions in two- 
dimensions, see Fig.[2j These will act as initial conditions 
for the coarsening dynamics. 



A. Area distributions 

A hull enclosed area is defined as the full interior of 
a domain boundary - irrespective of there being other 
interfaces and thus regions of the opposite phase within. 



The distribution of domain areas - we recall that do- 
mains are clusters of connected aligned spins - at critical 
percolation is given by Q 

n,iA, 0) ^ ^'"^3 , with r' = ^^ 2.055 . (2) 

Of course, the constants Cd and Aq are not separately 
defined by this relationship, only the ratio c^Ag In 
practice it is convenient to choose Cd to be the value ap- 
pearing in the domain area distribution at general times 
- see Eq. ([9]). The quantity Aq can be interpreted as a 
microscopic area, and will be introduced at various points 
as a small-area cut-off. The quantity Aq in Eq. (jSj) 
sets the units in such a way that [ud] = ^0 ^- This result 
is also valid in the limit Aa A <^ L^. 
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In equilibrium at Tc, Stella and Vanderzande [8| com- 
puted the number density of domains with area A 



A"- 



with r = 



379 
187 



2.027 (3) 



in the large A limit. Janke and Schakel [9] confirmed this 
claim numerically finding r « 2.0269. Motivated by the 
Cardy-Ziff result for hull enclosed areas, we conjecture 
that the prefactor Cd in ([3]) is the same Cd (up to terms 
of order c\) as that appearing in the prefactor 2cd for 
critical percolation. We discuss this point in detail at the 
beginning of section IV. B, and we check this hypothesis 
numerically in Sect. IVl 

We find it useful to include the small-area cut-ofl^, Aq, 
in these number densities, transforming the denominators 
to {A + Aq)"^ or {A + AoY''^ for hull enclosed and domain 
areas, respectively. 



Perimeter distributions 



equilibrium conditions using different analytic methods 
that include conformal invariance and renormalization 
group and Coulomb gas techniques 
Many numerical studies [1, [l^, [3, fTGl llTll confirmed and 
complemented the results in these analytic works. These 
works focused on the fractal dimensions of the domain 
area, Dd, and of the hull length, Dh- In two-dimensions 
these exponents are linked to the distribution exponents 
T and C/i as Q 



Dd 



r- 1 



Dh = 



a-1 



(6) 



and similarly for the primed quantities. 

We concentrate on the fractal properties of geometric 
structures by comparing the area of the clusters to their 
associated perimeter. This approach was used by Cam- 
bier and Nauenberg 18] who studied domain walls in the 
2dIM - with internal and external border - at equilibrium 
below, but near T^, and found 



One can equally study the length of hulls (external 
perimeters) and domain walls (including external and in- 
ternal boundaries). 

The number density nh{p,0) of hulls with length p, 
defined as the number of spins on the hull, at the critical 
point of the 2dlM was computed by Vanderzande and 
Stella [3 



A 



p 



1.43 



(7) 



"/i(p,0) 



Ch-3 



■ , . 27 
with C/i = Yi 



2.454 



(4) 



(P + Po)^'^ 

The value of the constant Cp^ was not estimated. As far 
as we know there is no analytic prediction for the num- 
ber density of domain walls at critical Ising conditions 
that should have the same functional form though with 
a possibly different exponent and a different constant 
■ Po is a microscopic length that we define as pp = Aq. 
For critical percolation, the number density of hulls 
with given perimeter was obtained by Saleur and Du- 
plantier (Tlj 



nh{p,0) 



(p +Po)^'' 



with — 



15 

y 



2.143 



(5) 



compatible with the numerical study |12l |. Again, the 
numerical value of the constant c' is not known. As far 
as we know, the analog of ^ for domain walls at critical 
percolation is not known, though we expect the same 
functional form with different constant Cp^ and exponent 

In Sect. IVII we show numerical results for the equilib- 
rium perimeter length number densities of critical Ising 
and infinite temperature configurations. 



C. Fractal properties 

Several authors studied the fractal properties of areas 
and perimeters in critical Ising and critical percolation 



The proportionality constant is not given. In Sect. lVII we 
revisit the geometric properties of the clusters at Tc and 
the corresponding one at very high temperature as well 
as their zero temperature evolution. 



III. SOME GENERAL RESULTS 

A number of general properties of hull enclosed and 
domain areas as well as hull and domain wall perimeters 
can be easily derived just by using the scaling hypothesis 
and two sum rules. We summarize them here. 



A. Scaling 

At long times, and irrespective of the initial condition, 
the total number of domain and hull enclosed areas per 
unit area, Nd.h{t) = /q dA nd.h{A,t), should scale as 
R~'^{t) in d dimensions, with R{t) a characteristic length 
scale usually associated to the 'typical' domain radius. 
For pure ferromagnetic coarsening, R{t) ~ (Xdi)^^"^ from 
which it follows that Nd^h(t) ^ in d = 2. Since char- 
acteristic areas scale as R^{t) ~ t, the scaling hypothesis 
implies that the domain and hull enclosed area distribu- 
tions have the forms nd.h{A,t) = t~'^ fd.h{A/t). 

In Sect. IV we present arguments that the initial dis- 
tributions of hull enclosed and domain areas determine 
the forms of these distributions at late times. In particu- 
lar, the initial forms are modified in a rather simple way 
at later times, such that the scaling forms are already 
suggested by the initial conditions. In the following we 
outline the consequences of this line of argument and de- 
fer detailed analysis to Sect. IV. 
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Let us first discuss domains with critical Ising initial 
conditions. Retaining the initial form ([3]), but includ- 
ing an extra time-dependent factor to satisfy the desired 
scaling form at large times, the expression for nd{A, t) is, 



nd{A,t) 



Cd [Xd{t + to)Y-'- 

[A + Xd{t + to)Y 



(8) 



at long times. Here, to is a microscopic time, and is 
a phenomenological parameter with the dimensions of a 
diffusion constant. Setting t = one recovers Eq. ([3|) 
provided Aq = Xdto and A^ Aq. Indeed, to is defined 
through this requirement while Aq is the by now usual 
microscopic area. For large t, Iq can be neglected and the 
conventional scaling form is recovered. 

For infinite temperature initial conditions we propose 



nd{A,t) 



2cd [Xd{t + to)]^'- 
[A + Xd{t + to)]-' 



(9) 



that also satisfies the scaling form asymptotically. 

One can easily check that the obvious generalization of 
Eq. |T]) including time 



nh{A, t) 



{A + Xntf 



(10) 



also has the desired scaling form. The constant Xh is 
another phenomenological parameter. Its value will turn 
out to be very close to A^. The factor 2 in the brackets 
in Eq. (fTO|) takes into account the two types of initial 
condition. 

We remark that in Eq. pUj) we do not explicitly include 
a short-time cut-off t^, while it is necessary in Eqs. ^ 
and (O to connect smoothly to the initial condition. 

Scaling can also be used to predict the time- 
dependence of the number density of hull and domain 
wall lengths. It yields 



nd,h{p, t) 



[p+{Xd,htY'^ +Po\'^-'-'^ 



(11) 



with the constant and exponent values depending on the 
initial condition and whether we are studying hulls or 
domain walls. This form is based on the assumption that 
the characteristic perimeters evolve in time as {Xdjit)^^"^- 



B. Sum rules 

We now present two exact sum rules which provide 
useful input for the analysis of nd,h{A^t). These sum 
rules apply at all times i > 0, for any initial condition, 
and for any working temperature. 

The first sum rule follows from the fact that the total 
domain area, per unit area of the system, is unity since 
each space point (or lattice site) belongs to one and only 
one domain. This gives 



dA A nd{A,t) = 1 . 



(12) 



The second sum rule follows from the fact that the 
total number of hull enclosed areas, Nh{t), is equal to 
the total number of domains, Nd{t), since each domain 
can be associated with a unique hull, namely the hull 
that forms its external boundary. This yields 



Nd{t) 



dAndiA,t) 
dAnh{A,t) = Nh{t) 



(13) 



Equations H])-® and their generalization describing 
the time-dependence, have been shown to hold for large 
areas (A 3> Aq) only and the number density can take a 
different form at small values of A. From the constraints 
(fn)l and (Uni), using Eqs. ©, ^ and we shall derive 
approximate relations between the constants c^, Ch-, Xd, 
Xfi, T and r' expected to hold at any working temperature 
T. These relations are exact to first order in the small 
quantity Ch- 

For the number densities of hulls or domain walls with 
given perimeter we have only one sum rule. The total 
number density must equal the number density of do- 
mains and hull enclosed areas: 



^Pd,ph (^) 



dp nd,h{p, t) = Nd{t) = Nh[t) . (14) 



1. Critical Ising initial conditions. 



The constraint on the total area using Eq. ([5]) for 
the number density of domain areas yields 

Cd = (t-2)(t- 1) « 0.02745 (15) 

where the numerical value was obtained for r = 379/187. 
If one takes into account that the minimal area is Aq (and 
not 0) the result is slightly different and it varies, though 
very weakly, with t. Indeed, the value of Cd decreases 
from Cd = {t - 2)(t - 1)2^-7^ at t = to Cd = (r - 
2)(t- 1) at t > <o = Ao/Xd- Inserting t = 379/187 gives 



0.02745 <cd< 0.02760, 



(16) 



a rather narrow interval. Note that Cd computed using 
P5p is quite close to the analytical result for c/i, namely 
c^" = l/87r\/3 ~ 0.02297. Indeed, Cd - c|" « 0.00448 
(which is order c^"/5) where we used the minimum value 
of Cd and the analytical value of Ch- 

We now consider the sum rule (|13p. at times t = Q 
and t 3> Ao/A/i, distinguishing the integral over the full 
interval [0,oo) from the one that takes into account the 
finite minimal area [Ao,oo). We denote by Ai the lower 
limit of the integration interval. We also include a small 
area cut-off Aq in the hull enclosed area number density, 
n,,(A,0)-c^/(A + Ao)2. 
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At time t = we find 



Cd 



Ch (r - 1) 



AT-2 
^0 



(17) 



Using now the expressions for Cd derived above we relate 
Ch to r: 



2) w 0.02674 
' 0.02639 



^1-0, 
Ai = Aq 



(18) 



These values are quite close to the analytic one c^" « 
0.02297. The relative difference is {ch - cr)/cr ~ 0.164 
{Ai = 0) and {ch - c^")/C « 0.149 {Ai = Aq). Note 
that using the sum rules we over-estimate the value of Ch- 
We can then expect to have over-estimated the value of 
Cd too. This remark will be important when comparing 
to numerical data. 

Finally, we can evaluate condition (|13p at late times, 
t :s> Aq/Xh, using the conjecture ([8]) for nd and the re- 
sult nh{A,t) ~ Ch/{A + XhtY deduced from scaling in 
Sect, nil XI and to be shown analytically in Sect. lIVXl In 
this case we find, independently of ^i: 



Ad _ Cd 1 
Xh Ch (r - 1) 



(19) 



This equation can be used to relate the two factors 
and Xh- Indeed, replacing Cd and Ch by their expressions 
as functions of r one finds 



Ad 
A,, 



1 , 

2^-2 



1.019 



Ai = 0, 
Ai = Ao 



(20) 



If, instead, we use Cd as derived above and the analytic 



0.022972 we find 



1.164 < < 1.170 

Xh 



(21) 



Since we derived Eq. (jl9p using the sum rules, it seems 
more appropriate to use the values of Cd and Ch obtained 
from the same relations. The values of Ad and Xh are then 
very close and consistent with the relation that we shall 
find in the next Section using the approximate equation 
for the evolution of domain areas [while Eq. pT|) yields a 
too large value for Ad]. 

The condition Np^{t) = Np^{t) = Nh{t) = Nd{t) im- 
plies 



Cp^ = -Ch (1 - C) , 
1-C 

Cpd — Cd , 
i — r 

for critical Ising initial conditions. 

2. Infinite-temperature initial conditions. 



(22) 



The infinite-temperature initial conditions turn out to 
be, after just a few time steps, equivalent to critical per- 
colation ones - see the numerical evidence in Sect. IV AI 
The sum rules yield, in this case. 



Cd 
Ch 

Ad 
Xh 



(r'-2)(r'-l) 



0.029 



Cd 



(r' - 1) 

Cd 

Ch (r' - 1) 



0.027 , 

= 1 , 



(23) 
(24) 

(25) 



where, for simplicity, we present results obtained with 
Ai = only. 

The conditions Np^{t) = Np^{t) = Nh{t) = Nd{t) im- 
ply 



Ph. 



= -2ch (1 - c;.) , 



- 2cd i-^ 



(26) 



for critical percolation. 



IV. STATISTICS OF AREAS: ANALYTIC 
RESULTS 

Our analytic results are obtained using a continuum 
description of domain growth in which the non conserved 
order parameter is a scalar field, (j){x,t), defined on a 
d-dimensional space. For a review in this problem, see 
[H. Its evolution is determined by the time-dependent 
Ginzburg-Landau equation or model A dynamics: 



7- 



d(j){x, t) 
dt 



svi<t>) 

S<p{x, t) 



+ e(f,t). (27) 



The potential y is a symmetric double well, with V{4) — > 
±oo) = oo and two minima at ±(/)o. ^ is a Gaussian 
distributed random scalar field with zero mean and cor- 
relation 

{i{x,t)£,{x' ,t')) = 2k^T-f5'^{x-x')5{t-t') . (28) 

This white noise introduces thermal agitation. T is the 
temperature of the thermal bath, is the Boltzmann 
constant and 7 is the friction coefficient. From now on 
we set the units in such a way that fee = 7 = 1 . The low- 
temperature ordering dynamics from a disordered initial 
condition corresponds to the growth of ordered domains 
of the two equilibrium states, (/)(x, t) = ±00, separated by 
interfaces. Using the evolution equation (j27p at zero tem- 
perature, Allen and Gahn showed that in any dimension 
d the velocity, w, of each element of a domain boundary 
is proportional to the local interfacial mean curvature. 



Xh 
'2tt 



(29) 



Xh is a material constant with the dimensions of a dif- 
fusion constant, and the factor l/27r is for later conve- 
nience. The velocity is normal to the interface and points 
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in the direction of reducing the curvature. The dynamics 
is then purely curvature driven at zero temperature. 

Temperature fluctuations have a two- fold effect. On 
the one hand they generate equilibrium thermal domains 
that are not related to the coarsening process. On the 
other hand they roughen the domain walls thus opposing 
the curvature driven growth and slowing it down. Then 
equation is no longer valid. However, it has been 
conjectured and verified numerically that, at least for av- 
eraged dynamic quantities well-described by the scaling 
hypothesis, all temperature effects are captured by intro- 
ducing a T dependent parameter. We shall use this 
working hypothesis in the analytic part of our paper and 
we shall put it to the test numerically. 

The number density of hull enclosed or domain areas 
at time t as a function of their initial distribution is 



n{A,t)= dA,S{A-A{t,A,))n{A„U) 



(30) 



with Ai the initial area and n{Ai,ti) their number dis- 
tribution at the initial time U. A{t,Ai) is the hull en- 
closed/domain area, at time t, having started from an 
area Ai at time U. 



A. Hull enclosed areas 

In two dimensions [20*1 we can immediately deduce the 
time-dependence of the area contained within any finite 
hull by integrating the velocity around the hull: 



dA 



V al = 

27r 



hidl ~ —Xh 



(31) 



In the second equality we used the zero-temperature 
Allen-Cahn Eq. and in the final one we used the 

Gauss-Bonnet theorem. Integrating over time, with ini- 
tial time tj, we find 



Ah{t, A,) = A., - Xh{t - U) 



(32) 



Therefore 

nh{A,t) 



dA^ S{A - A^ + Xh{t - U)) nh{A,,ti) 



nh{A + Xhit-U),ti) 



(33) 



In deriving this result we have implicitly assumed that 
a single domain cannot split into two, and that two do- 
mains cannot coalesce. A little thought shows that nei- 
ther process is possible for two-dimensional curvature- 
driven growth since both processes require that two parts 
of a single domain boundary (for splitting) or parts of 
two different domain boundaries (for coalescence) come 
together and touch. But it is clear that the curvature 
driven dynamics always acts to prevent this happening, 
since the velocities of the domain boundaries at the in- 
cipient contact point are in opposite directions. 



The initial distributions, nh{Ai,ti), are given by the 
Cardy-Ziff results displayed in Eq. ([T|) - assuming, in 
the case of a quench from infinite temperature, that the 
system rapidly sets into the critical percolation condition 
(see Sect. IV A| . For t:s> ti one immediately recovers the 
results in 



nhiA,t) = 
nh{A,t) = 



2ch 



[A + Xhtr 

Ch 

{A + Xhtf 



To ^ oo , (34) 
To = T, , (35) 



in the limit Aq <^ A <^ , i.e. for hull enclosed ar- 
eas much larger than microscopic areas but much smaller 
than the area of the system. 

Equations and ([55]) have the expected scaling 

forms nh{A,t) ^ t^'^f{A/t) corresponding to a system 
with characteristic area proportional to t or characteris- 
tic length scale R{t) (x t^^^, which is the known result 
if scaling is assumed Here, however, we do not as- 
sume scaling - rather, it emerges from the calculation. 
Furthermore, the conventional scaling phenomenology is 
restricted to the 'scaling limit': A — > oo, t oo with 
A/t fixed. Equations ([M]) and ([55]) . by contrast, are valid 
whenever t is sufficiently large and A^ Aq. This follows 
from the fact that, for large t, the forms and ([551) 
probe, for any A^ Aq, the tail (i.e. the large-A regime) 
of the Cardy-Ziff results, which is just the regime in which 
the latter is valid. The restriction A ^ Aq is needed to 
justify the use of Eq. which breaks down when the 
reciprocal of the curvature becomes comparable with the 
width of a domain wall. 

The averaged area enclosed by a hull is then given by 



_ fdA'A'n,iA\t) 
cx Xht 



(36) 
(37) 



with a time-independent prefactor that behaves as 
(AqIuL^) for large system sizes. The reason for the di- 
vergent prefactor in the infinite size limit is that a site 
can belong to several hulls. 



B. Domains 

For the domains we need to write an evolution equation 
and derive, at least approximately, the area at time t, 
Ad{t,Ai), of a domain with initial area Ai. We shall 
show that the time-dependent number density of domain 
areas is indeed given by our guess, Eqs. ([8]) and ([9]), for 
the two classes of initial conditions. 

Our strategy is to exploit the smallness of the parame- 
ter Ch « 0.023. Although Ch is a constant, we can exploit 
a formal expansion in Ch in the following sense. Since the 
total number of hulls per unit area is proportional to Ch, 
the number of interior hulls within a given hull is also pro- 
portional to Ch, and so on. This means that, in dealing 
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with domains we need consider only the first generation 
of interior huUs, since the number of "hulls within hulls" 
is smaller by a factor Ch- With this approach, only one 
approximation - a kind of mean-field one on the number 
of first-generation hulls within a parent hull (see below) 
- is necessary. 

The same line of reasoning shows that, in a hypothet- 
ical theory in which Ch can be treated as variable, the 
distinction between hulls and domains will disappear in 
the limit Ch — *■ 0. In this limit, therefore, the exponents 
T and r' must both approach the value 2, i.e. we can 
formally write r = 2 + 0{ch) and r' = 2 + 0{ch)- Fur- 
thermore, due to the factor 2 that appears in but 
not in ((35| . the ratio (r' — 2)/(t — 2) must approach the 
value 2 in the limit — > 0. The actual value of this ratio 
is 187/91 = t' = 2.055, not very far from 2. Indeed the 
difference is of order Ch as expected. 

We can use the same line of argument to discuss c^, Ch, 
Xd and Xh- Since in the (hypothetical) limit Ch — > 0, hulls 
and domains become identical, it follows that in this limit 
one must have Cd Ch and Xd Xh, i.e. Cd = Ch + 0{cD, 
and Xd — Xh[l + 0{ch)]- AH of these results are consistent 
with the relations (fT5|) and (fT9|) derived from the sum 
rules ((El) and (dSl). 



The evolution of domain areas. 



Take a hull with enclosed area Ah at time t. This hull 
is also the external border of a domain, which may it- 
self contain one or more 'first level' sub-domains whose 
external borders form the internal border (which may 
be disconnected) of the original domain. These external 
borders of the first level sub-domains are themselves 'first 
generation' hulls lying within the parent hull. These in- 
terior hulls can themselves have interfaces in their bulk 
separating domains of the reversed phase (higher gener- 
ation hulls), see Fig. [3] where we show a sketch with this 
structure. 




FIG. 3: (Colour online.) Sketch of a configuration with four 
circular hulls and domains. The parent hull has radius _Ri. 
There are two first generation hulls with radius R2 and Rs 
and one second generation hull with radius -R4. = 2 in 
this example. The interior border of the external domain is 
disconnected and has two components. 



Let us call i'{t) the number of first-generation hulls 
within the parent one. It is clear that i>{t) is semi-positive 
definite, monotonically decreasing as a function of time 
and reaching zero at a given instant tmax , when all inte- 
rior hulls disappear and Ad — Ah thereafter. One can es- 
timate tmax from = A''j^\tynax) = ^^'l^O ^ ^'hitraax - 

ti), which yields t^ax — ti= ^h'''^ (ti) / Xh where the index 
int indicates that we are studying here the first gener- 
ation hull with maximal initial area (all others having 
already disappeared) . It is clear that tmax — U is smaller 
but of the order of Ah{ti)/Xh where we replaced A™*(ti) 
by the initial area of the parent hull: 



(t. 



AhjU) 
Xh 



(38) 



We wish to write a differential equation for the time- 
evolution of the parent domain area. It is clear that, at 
first order in dt: 



Ad{t + dt) = Ad{t) - Xhdt + iy{t)Xhdt 



(39) 



where the second term in the right-hand-side represents 
the loss in area due to the inward motion of the exter- 
nal domain-wall while the last term is the gain in area 
due to the outward motion of the first-generation internal 
domain- walls. This gives 



dAdjt) 
dt 



-Xh [1 - Ht)] 



(40) 



Differently from hull enclosed areas, that always de- 
crease in size as time passes, domain can either diminish 
{i^ — 0), increase {ly > 1) or conserve {i^ = 1) their area 
in time. 



2. The number of first generation interior hulls. 



We cannot, of course, know the exact number of first 
generation hulls falling within a selected hull with en- 
closed area Ah- We can, however, estimate it with an 
upper bound obtained by counting all interior hulls and 
averaging over all parent hulls using nh{A,t) derived in 
Sect. IIV Thus, we expect 

iy{t) < {v{t))A^(t) 
rAh(t) 

{^{i))A,At) - Mt) / dAnh{A,t) 



ChAl{t)[Xh{t - t,) + Aq]- 
[Ah{t) + Xh{t - t^) + Ao] 



(41) 



where we include a small area cut-off, Aq , in the denomi- 
nator of Uh and, for concreteness, we use the hull enclosed 
area distribution for critical Ising initial conditions. This 
equation can be further simplified if one uses that at time 
t the hull enclosed area we are interested in is given by 



Ahit) ^ Ahit^) ~ Xhit - t,) 



(42) 
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(We call here Ah{ti) the initial area of the hull.) Then 
Ch[Ah{U) - Xhit - U)]'' 



'At) 



[Xh{t-U) + Ao] [AhiU) + Ao] 



(43) 



Note that, although we over-counted the interior hulls by 
including second-generation, third-generation, etc. hulls, 
the number of these is of order c^, c^, . . . respectively, 
so this treatment is exact to leading order in Ch except 
for the replacement of i^(t) by its average over all first- 
generation hulls of the same area. 

The most interesting cases are such that Ah{ti) ^ Aq, 
otherwise the hull and domain areas are just identical or 
very similar. In these cases ( iy{ti) )Ahiti) ^ ChAh{ti) / Aq. 
Expression (I43p has the following limiting values 



>(*)>A,(t) 



Ch Ah{ti) 
\h{t-U) + Aa' 



ach, 



AhiU) > Xhit-U) 



Ah{U) - Xh{t-U) 



we used Ah{ti) S> Aq in the last case, and a is a numerical 
constant of the order of Ah{ti). The result is a very 
small quantity, of the order of Ch, in both cases. The 
remaining mathematical possibility, Ah{ti) < Xh{t — ti) 
is not realized because Ah{t) cannot be negative. 

While u{t) vanishes at tmax^ see Eq. (1551) . ( i'{t) )Ah(t) is 
different from zero at all times. Thus, Eq. (|43)) cannot be 
used beyond the limit t„iax when all internal hulls have 
already disappeared and it is no longer correct to replace 
v{t) by {v{t))A,[t)- 

The analysis of infinite temperature initial conditions 
is identical to the one above with Ch replaced by 2ch- 
Thus, {i^{t)) Ahit) is expected to take twice the value it 
takes for critical Ising initial configurations. 

We have checked the accuracy of this approximation 
numerically by counting the number of first generation 
internal hulls falling within each parent hull at different 
times. Figure m shows the results for the zero tempera- 
ture evolution of the 2dIM starting from Tq — > cxd and 
To = Tc initial conditions. While at very short times one 
sees deviations between the numerical data and analytic 
prediction, the agreement between the two becomes very 
satisfactory for times of the order of i = 64 MCs and 
longer, as shown in the figure. 



The instantaneous domain area. 



If we now replace h'{t) by ( h'{t) )AH(t) given in Eq. (|^ . 
it is quite simple to integrate the differential equa- 
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FIG. 4: (Colour online.) Comparison between p{t) and 
{ u{t) )A^(t) for the r = evolution of the 2dIM with To ^ oo 
and To = Tc initial conditions. The measuring time is 
t — QA MCs. The curves are given by Eq. (|4T} in the 
limit t ^ ti and Ah(t) 3> ^o, leading to the functional 
form {i'{t)) A^(t) = Ch2^^/(1 + s) ~ ChX when x ^ 1, with 
X = Ah{t)/Xht and Xh = 2.1, see Sect.|V]for an explanation 
of the choice of this value. 



tion (Sol). One finds 



1 2 



Adit) = AdiU) - Xhil + 2ch)it ~ t,) 



ChAl 



2[Ah{U)+AQ] 
- Ch[Ah{ti) + Ao\\tl 



1 



\h{t-ti) 



Ao 
Xhjt-U) 
Ao 



Setting t = ti one recovers Ad{t) — Ad{ti) as required. 
In the natural cases in which Ah(ti) ^ and for long 
times such that Xh{t — ti) 3> Aq this expression can be 
rewritten as 



Ad{t) = Ad{ti)-Xh 
+ ChAh{ti)\n 



1 



Ch Kit ~ ti) 



2 Ahiti 
Xhit-t,)' 



Ao 



The factor in the second term 



Xdit) = Xh 



l + 2ch 



Ch Xhit - ti) 



2 Ahit,) 



it - ti) 
(44) 

(45) 



is a very weakly time-dependent function. Since t can 
take values between the initial time, t = ti, and the max- 
imum time before the first generation hull itself disap- 
pears, tmax = ti + Ahiti) /Xh, Xdit) varies within the 
interval: 



A. 1 



3ch 
2 



< Xdit) < Xhil + 2ch) 



(46) 



These bounds are indeed very close. As expected from 
the analysis of the sum rules, see Sect. Illli takes a 
slightly higher value than A^; it equals A;i plus a small 
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correction of order cu (in practice, 1.035 A/i < < 
1.046 A/i using the analytic value for Ch)- 

The coefficient in front of the logarithm, ChAh{ti), is 
0{cfi). The sum rules, imply Ch = Cd + 0{c\). Neglecting 
the higher order correction we can then replace Ch by Cd- 
The same applies to Ah{ti), which equals Ad(ti) plus a 
term 0{ch) that we can equally neglect. Thus Ah{ti) ~ 
Ad[ti) = Ai. In this way we obtain 



Ad{t, A,)~ Ai-\d[t- U) + CdAi In 



1 



Aq 



Inserting this result into Eq. (|30p . including the micro- 
scopic area Aq as a small-area cut-off in the denominator 
of Eq. ID) , then gives 



nd{A,t) ~ CdAl 2 J 1 + Cdln 



A, 



T-l 



x[A + Xd{t-t,) + Ao]-^ 



(47) 



where we have replaced Xh by A^ inside the logarithm, 
which is correct to leading order in Ch- Using the fact 
that Cd is very small and of the order of (t — 2)(t — 1) = 
(r — 2) +0{cf^), as implied by the sum rules, we can now 
exponentiate, correct to leading order in c^, the logarithm 
in the curly brackets to obtain 



ndiA,t) 



Cd [Aq + Xdjt - t,)]-' 
[A + Ao + Xdit - U 



(48) 



Finally we set the initial time, to zero and write the 
microscopic area, Aq, as Xdto to obtain the expected form 



nd{A,t) 



Cd [Xdjt + tQ)Y-'- 

[A + Xd{t + ta)Y 



(49) 



for the time-dependent number density of domain areas. 

The same sequence of steps for infinite-temperature 
initial conditions leads to the same form but with Cd re- 
placed by 2cd and r replaced by r'. The effects of tem- 
perature are expected to appear only through the param- 
eters Xd and A/i once thermal fluctuations are extracted 
from the analysis. 

The averaged domain area is then given by 



{A){t) 



1 



oc 



J dA'A'ndiA',t) 

f dA'nd{A',t) " N^tj 
Xdt . 



(50) 
(51) 



V. 



STATISTICS OF AREAS: 
TESTS 



NUMERICAL 



To test our analytic results we carried out numerical 
simulations on the 2d square-lattice Ising model {2dlM) 
with periodic boundary conditions using a heat-bath al- 
gorithm with random sequential updates. All data have 



been obtained using systems with size = 10'^ x 10'^ 
and 2 x lO'^ runs using independent initial conditions. 

Domain areas are identified with the Hoshcn- 
Kopelman algorithm ^2T| while huU-cncloscd ones are 
measured by performing a directed walk along the inter- 
faces, in analogy with the algorithm in [3]. A detailed 
description of our algorithm is given in Appendix [XI 

The equilibrium critical Ising initial conditions have 
one spanning cluster (since the system is at the limit of 
the percolation threshold), that grows during evolution. 
No other spanning cluster is later formed. By contrast, 
equilibrium infinite temperature initial conditions are be- 
low the critical random percolation point in d = 2 but 
often after 2 MC steps two spanning clusters appear that 
also grow during evolution. After 20 MCs roughly 50% 
of the spins lie typically on these clusters. This implies 
that we need to simulate a large number of independent 
samples to obtain a good statistics. 

It is important to note that the dynamics of the dis- 
crete model includes processes that are not taken into ac- 
count in the continuous model, as given in Eq. (|29p . for 
which we derived our analytical results. Some of these 
processes are the fission of a big domain into two smaller 
ones (that usually occurs by cutting a thin bottle neck 
that joined them), or the coalescence of two domains to 
form a bigger one. However, we shall prove that these 
processes are not important and the dynamics of the dis- 
crete model is well described by the analytic results. 



A. Initial conditions. 

We used three types of initial conditions: equilibrium 
at infinite temperature, Tq — s- cx); equilibrium at the crit- 
ical point, Tq = Tc] equilibrium within the high temper- 
ature phase, To = Tc + AT with AT > 0. 

We mimicked an instantaneous quench from Tq ^ oo 
by using random initial states with spins pointing up or 
down with probability 1/2. Assigning site occupation to 
up spins and vacant sites to down spins the infinite tem- 
perature initial condition can be interpreted as a per- 
colation problem a,t p = 0.5 and thus below the ran- 
dom percolation transition pc = 0.5927 in a square bi- 
dimensional lattice. Even if initially away from critical- 
ity, in a few MC steps the hull enclosed area distribution 
becomes the one in Eq. ([1]), as shown in the top panel 
in Fig. [5l The initial distribution lacks large areas, there 
being almost none with A > 10^, and the tail of Uh falls 
off too quickly well below the critical percolation curve. 
In a few time steps large structures appear and the tail of 
the distribution approaches the expected form at critical 
percolation. Simultaneously, the weight at small areas 
diminishes and the curve progressively gets flatter. This 
effect can also be seen in Fig. [55] where we display data 
for perimeter lengths. In the bottom panel in Fig. [5| we 
display our numerical results for the domain area distri- 
bution, which are compatible with the form It is 
intuitively clear why this must be so. If the system is 
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A 

FIG. 5; (Colour online.) Early evolution of the infinite tem- 
perature initial condition. Top panel: in a few MC steps the 
hull enclosed area distribution reaches the one of critical per- 
colation lU, whose slope, -2, is shown by the straight line 0|. 
Bottom panel: the domain area disjtribution. The straight 
line is the power-law decay ^-2.055 



coarse-gramed on the domain typical scale, R{t), it will 
look completely disordered. When R{t) is large compared 
to the lattice spacing, the disorder will be that of contin- 
uum percolation, for which the critical density is one half 
by symmetry in two dimensions [s^ . [33| . It follows that 
the coarsening system will be asymptotically at percola- 
tive criticality, i.e. the dynamics self-tunes the system to 
percolative criticality in two dimensions [provided R{t) 
remains much smaller than the system size]. The data 
show that, as far as the hull and domain area distribu- 
tions are concerned, this only takes a few Monte Carlo 
steps in practice. During these few steps many small 
domains coalesce to form larger ones meaning that the 
dynamics is dominated by processes that are not taken 
into account by Eq. ((29)l . This argument also shows that 
the domain distribution in the scaling limit indeed has 
the predicted A^'^ tail. It is interesting to remark that 
the system approaches the percolative critical state not 
by increasing p (indeed, the magnetization during the 
coarsening process initially remains close to zero), but 
by decreasing the value of Pc, from 0.5927 to 0.5, as the 
correlation between spins increases. 



We can look at this from another perspective in the 
context of the continuum model. Consider a random 
field (/)(x), symetrically distributed with respect to zero, 
with bounded variance and two-point correlator C(r) = 
( <j){x)(j){x -f r) ) with r = |r^. The zero contour lines of 
this field can be imagined to divide the plane into re- 
gions of black and white with each contour line forming 
a boundary between black and white regions. Provided 
that C(r) falls off faster than r^^/^ for large r, this prob- 
lem is known to belong to the standard percolation uni- 
versality class [3^ . If we now identify 4>{x) with the order 
parameter field when well-defined domain walls (the zero 
contour lines) have formed, we see that the resulting do- 
main structure corresponds to critical percolation. 

In the plots we use a double logarithmic scale that 
serves as a first check of the power-law decay of the prob- 
ability distributions but it is not accurate enough to ex- 
amine the value of the constants cu and Cd- We delay the 
presentation of a very precise test of these parameters 
to Sect. IV Bl where we analyse the time-evolution of the 
distribution functions. 

We obtained the initial states for the coarsening dy- 
namics at the critical temperature, Tq = Tc, and at 
To = Tc AT, after running 10^ Swendsen-Wang al- 
gorithm steps. We checked that the systems are well 
equilibrated after these runs. The distribution of hull 
enclosed and domain areas at Tc are consistent with the 
analytic forms ~ not shown. 




A 



FIG. 6: (Colour online.) Equilibrium number density of hulls 
per unit area for the high temperature phase of the 2dIM. The 
distributions seem power law for small areas, with a temper- 
ature dependent exponent that approaches 2 when T — > Tc. 
Indeed, data for T = 2.3 ~ 1.01 Tc is almost coincident with 
Eq. JSSl). 
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B. Coarsening at zero temperature 

1. Hull enclosed areas. 

In Fig. [7] we show the time-dependent huh enclosed 
area distribution in double logarithmic scale, at seven 
different times, following a quench from Tq ~* oo- The 
figure shows a strong time dependence at small areas and 
a very weak one on the tail, which is clearly very close to 
a power law. The curves at small areas move downwards 
and the breaking point from the asymptotic power law 
decay moves towards larger values of A for increasing t. 




A 



FIG. 7: (Colour online.) Number density of hull enclosed ar- 
eas per unit system area for the zero-temperature dynamics 
of the 2dIM at seven times evolving from an infinite tem- 
perature initial condition. The lines represent Eq. (|34p with 
Ch = l/87r^/3. 
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FIG. 8: (Colour online.) Number density of hull enclosed ar- 
eas per unit system area for the zero temperature dynamics 
of the 2dIM evolving from an infinite temperature initial con- 
dition. The full line is the prediction (|34p with = l/87r\/3 
and Ah = 2.1. 

In Fig. [TO] (bottom) we zoom on the small area region 
{A < 10^) where the time-dependence is clearer and we 
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FIG. 9: (Colour online.) Finite size effects at t = 16 MCs; four 
hnear sizes of the sample are used and indicated by the data- 
points. The value of A/ 1 at which the data separate from the 
master curve grows very fast with L with an exponent close 
to 2. 



scale the data by plotting {Xiit)'^nh{A,t) against A/Xht 
with A;i — 2.1. We tried other time-dependent factors but 
Xht with this particular value of Xh is the one yielding the 
best collapse of data at small areas, Aq ^ A < Xht. For 
A larger than the 'typical' value Xht the time and de- 
pendence becomes less and less important. In Fig. |S] we 
show the data in their full range of variation in log-log 
form to test the prediction nh{A,t) oc A~'^ for large A. 
The data are in remarkably good agreement with the pre- 
diction p4p ~ shown as a continuous curve in the figure 
- over the whole range of A and t. The downward devia- 
tions from the scaling curve are due to finite-size effects. 
The latter are shown in more detail in Fig. [SI where we 
display the t = 16 MCs results for several linear sizes. 
Finite size effects appear only when the weight of the 
distribution has fallen by many orders of magnitude (7 
for a system with L — 10^) and are thus quite irrele- 
vant. In the tail of the probability distribution function 
(pdf) the numerical error is smaller than the size of the 
data points. The nearly perfect agreement between the 
analytical theory and the data is all the more impres- 
sive given that the curvature-driven growth underlying 
the prediction ((M)) only holds in a statistical sense for 
the lattice Ising model [22|. Indeed, even at small values 
of A/Xht, where the lattice and continuous descriptions 
are expected to differ most, the difference is only a few 
percent, as we shall show below. 

It is clear that the evolution of the hull-enclosed area 
distribution follows the same 'advection law' (j33p . with 
the same value of Xh, for other initial conditions. The 
evolution from a critical temperature initial condition is 
shown in the top panel of Fig. [TOl A fit of the data at 
small areas yields the value of the parameter Xh that, 
consistently with the analytic prediction, takes the same 
value A/i = 2.1. In the bottom panel, we compare the 
time-dependent hull enclosed area distributions for the 
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initial conditions Tq — > cx) and Tq = Tc and we zoom 
on the behaviour of nh{A,t) at small areas, A/XyJ: < 10. 
The two solid lines correspond to the numerator in 
being equal to 2ch for infinite temperature initial con- 
ditions and c/i for critical Ising initial conditions. The 
difference between the numerical data for the two initial 
states is clear and it goes in the direction of the analytic 
prediction (a factor 2 difference in the constant). Finally, 
while the log-log plot in Fig. [8] suggests that the data are 
compatible with ~ 0.023 this way of presenting the 
data is not precise enough to let us quantify the accu- 
racy with which we match the analytic prediction. We 
test the numerical values of the constant Ch in detail in 
Sect. IVl where the numerical error is also estimated. 
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conditions at different Tq > Tc show only a different tran- 
sient behaviour: the closer they are from Tc, the longer 
it takes to reach the asymptotic law, Eq. Equilib- 
rium initial distributions, for Tc < Tq < oo, are shown 
in Fig. [6l while Fig. 1111 shows an example of their subse- 
quent evolution. In the latter, both analytic predictions, 
for To = Tc and Tq = oo, are shown as solid lines along 
with data for increasing times after a zero-temperature 
quench from Tq = 2.5. In the first steps, the curve fol- 
lows the one for critical initial conditions at small A/Xht 
and then departs to reach the one for infinite tempera- 
ture initial conditions at large A/Xht. At longer times, 
the deviation from the critical initial condition line oc- 
curs at a smaller value of A/Xht. Initially the system has 
a finite, though relatively small, correlation length £,{T). 
Thermal fluctuations with linear size of the order of ^ 
and also significantly larger than ^ exist (see the discus- 
sion on the effect of thermal fluctuations in Sect. IV Cp . 
Notice that £,{T) does not correspond exactly to the size 
of geometric domains: thermal fluctuations are not per- 
fectly described by domains of aligned spins, since not all 
of them are correlated. At any given temperature above 
Tc, fluctuations smaller than £,{T) have the same statis- 
tics than those occurring at Tc and are thus described by 
Eq. ([55)) . while domains larger than £,{T) are not made 
of correlated spins and thus are described by the infinite 
temperature distribution, Eq. (j34p . As time increases, 
the system loses memory of the finite-size fluctuations 
and the asymptotic state does not differ from the infinite 
temperature one. Only when fiuctuations exist over all 
spatial scales does the asymptotic state differ. 

This behaviour can be interpreted as follows. At fixed 
A/t, shorter times correspond to small areas while longer 
times are related to larger areas. Very small areas corre- 
spond to short linear sizes, of the order of the domains 
in the initial configurations, and thus reminiscent of crit- 
ical ones. Instead large areas correspond to long linear 
sizes that are much longer than the correlation length and 
closer to the ones reached from the infinite temperature 
initial condition. 



2. Domain areas. 



FIG. 10: (Colour online.) Number density of hulls per unit 
area for the zero-temperature 2dIM evolving from critical ini- 
tial conditions. The initial states are obtained after running 
10^ Swendsen-Wang algorithm steps. Top panel: the full (red) 
line is Eq. (|35|l with Xh ~ 2.1 which again yields the best fit 
of the data at small areas. For comparison we include with a 
dotted (blue) line the analytic prediction for an infinite tem- 
perature initial condition, i.e. Eq. ((34}. Bottom panel: more 
details on the influence of the initial conditions. The two data 
sets correspond to configurations taken at several times after 
a quench from To — > oo and To = Tc. The solid lines are the 
analytic predictions (|34|l . blue line, and ([35}, red line. 

Moreover, Eq. ([M)l applies to any Tq > Tc equilib- 
rium initial condition asymptotically. Equilibrium initial 



We now study the domain areas and perform the same 
tests as above though focusing on the analytic predictions 
® and dH). 

In Fig.[T3we display the number density of domain ar- 
eas in the scaled form for two initial conditions, Tq — > oo 
and To = Tc, after removing any spanning domain from 
the statistics. For comparison, in Fig. [131 the same distri- 
butions with the spanning domains are shown. As done 
for the hull enclosed areas we fit the parameter Xd by 
analysing the behaviour at small areas, A < Xdt, and we 
find, once again that — 2.1 yields the best collapse of 
data (see the discussion in Sect. lIVBp . We use the ex- 
trapolated value Cd = 0.025 obtained with the numerical 
analysis described in full detail in Sect. IV 531 Note that 
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FIG. 11: (Colour online.) Effect of the initial condition tem- 
perature on the hull enclosed area distribution. The contin- 
uous lines are the analytic results for equilibrium initial con- 
ditions at To = oo and To = Tc, bottom (red) and top (blue) 
lines, respectively. In between, we present numerical data for 
two different times given in the key after the quench from an 
initial state equilibrated at To = 2.5. 



we expect the difference between ca and Ch to be of order 
c^, and thus rather hard to observe numerically. 

Both sets of figures, [12] and [131 exhibit finite size ef- 
fects in the tail of the distributions, where the number 
of domain areas has already decreased by several orders 
of magnitude. As for the hull enclosed areas, the point 
where these finite size effects cause the deviation from 
the collapsed curve moves towards the right as the sys- 
tem size increases, becoming less and less relevant. In 
Fig. [121 large domain areas (violating the limiting con- 
dition A<^1?^ that would nonetheless be accounted for 
in an infinite system are here removed since they span 
the system in one of the directions, leading to the down- 
ward bending of the distribution. On the other hand, in 
Fig. [131 when counting these domains, they are chopped 
by the system boundaries, thus contributing to the dis- 
tribution in a region shifted to the left, accounting for 
the bumps seen in the figure. 



3. Study of the constants and Cd . 



In order to improve the data analysis we followed the 
procedure used by Cardy and Ziff, who studied the 'fi- 
nite area' scaling of the cumulative distribution between 
A and 2 A [6]. The method is as follows. For hull enclosed 
and domain areas the total number of areas between A 
and 2A is Nh,d{A,2A) = f/ dA' nhA^'). Using the 
analytic prediction for nu one finds 2ANh{A,2A) ~ ct 
for large areas at critical Ising conditions and without 
the factor two at critical percolation. Following Ref. @, 
we assume that there are power-law finite area correc- 
tions and add a term like aA^^ to the above expressions. 
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FIG. 12: (Colour online.) Number density of domains per unit 
area for the zero-temperature 2dIM evolving from To oo 
(top) and To — Tc (bottom) initial conditions. In both figures 
the spanning clusters have been extracted from the analysis 
(compare with Fig. [TJIwhere we include them) . The full (red) 
line represents Eq. (|49ll . with Cd = 0.025 and r' = 2.055, 
Cd Cd/2 and t = 2.027 (bottom panel), and \d = 2.1 in 
both cases. 



From this relation one extracts the value of c/j. Similarly, 
for domains one can use (1 - T)Nd{A,2A)/[{2A)'^^'' - 
A^^^] = Cd + aA~^ at critical Ising initial and its modi- 
fied form at critical percolation initial conditions. 
In the dynamic case we have 



2A 

dA' nh,d{A',t) . 



Using the predictions for n^.d we find 

2chA 



Nh{A,2A-t) = 



and 



Nd{A, 2A- 1) 



{A + \ht){2A + Xht) 
'^^^'^'^"'\i2A^Xdt) 



(52) 



(53) 



l-r' 



1 -r' 

-{A + Xaty--'] (54) 



for Tq ^ cxD and without the factor two and with the 
exponent r' replaced by t for Tq — Tc- To extract the 
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FIG. 13: (Colour online.) The same as Fig. [TH but with the 
spanning domains included in the statistics. Notice that due 
to the finiteness of the lattice, the actual area of these span- 
ning domains are smaller than the value they would have on 
an infinite system, generating the overshoot of the distribu- 
tion for values of A close to . The larger is the system, the 
more to the right these peaks are localized. 



values of the constants ct and c^, we rewrite these forms 
as 



{2A)-\A + \ht){2A + Xht)Nh{A, 2A- 1) 



Ch + a 



A 



-6 



(55) 



and 



{Xdtf-^'{l-T')Nd{A2A;t) 
2[{2A + Xdty-^' - (A + Xdty-^ 



= Cd 



A_ 

Xd 



(56) 



and similarly for Tq — Tc- a is a constant that takes 
different values for different times. In Figs. fMlfTHl we show 
the outcome of this analysis. We use the same scale on 
the vertical axis in all plots to compare the accuracy of 
the results. 

In Fig. [T3] we show results for the hull enclosed areas. 
In the top panel we study equilibrium data at Tc- The 
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FIG. 14: (Colour online.) Study of Ch using finite area scal- 
ing implemented as in The numerical data correspond to 
2ANh{A,2A) in the equilibrium case and the left-hand-side 
of Eq. (|56|l in the dynamic case. Top panel: equilibrium dis- 
tribution of hull-enclosed areas at criticality. The upper hori- 
zontal line is the predictions for Ch stemming from the use of 
the sum rule. The lower horizontal line is the analytic predic- 
tion. The intercept of the inclined straight line with a; = is 
the numerical test of the analytic prediction. Bottom panel: 
study of Ch during coarsening after a quench from criticality 
to zero temperature. Note that the prediction from the fit is 
slightly higher than the analytic value, see the discussion in 
the text. 



extrapolation of the numerical data approaches the ana- 
lytic prediction for Ch within a 2% accuracy. We use the 
value of the exponent b proposed by Cardy and Ziff Q. 
The estimation from the approximate use of the sum 
rules is well above the analytic prediction and numeri- 
cal value. In the bottom panel we test the value of Ch 
in our dynamic prediction by studying data at three in- 
stants, t = 16, 64, 128 MCs, after a critical Ising initial 
configuration. In this analysis we use the same exponent, 
h = 0.875, as in the study of the equilibrium data and 
we fit the other parameters, Ch and a. The fits of the 
time-dependent data to straight lines extrapolate to the 
same value that is, however, of the order of 5% off the 
analytic prediction. Note the non-monotonic character of 
the time-dependence in the slopes a. One has to keep in 
mind that the extrapolated value is very sensitive to the 
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FIG. 15: (Colour online.) Study of during coarsening after 
a quench from To — > oo to T = using finite area scaling 
implemented as in [^. Note that the prediction from the fit 
weakly depends on time. 



fit, in particular, to how many data points are considered. 

In the analysis of the infinite temperature initial con- 
ditions we are forced to use dynamic results to reach, 
first the critical percolation situation, and next follow 
the coarsening evolution. In Fig.[Tn]we display this type 
of data for the instants given in the key. The prediction 
from the fit is slightly different from the analytic result at 
short times but it approaches the analytic value, shown 
with a horizontal line, at sufficiently long times (look at 
the t — 128 MCs results). In conclusion we find 



Ch ^ 0.0229 ± 0.0015 . 



(57) 
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FIG. 16: (Colour online.) Study of Cd using finite area scal- 
ing implemented as in [^. Cd extracted from time-dependent 
pdfs evolving at T = after a quench from criticality. The 
numerical value of Cd is roughly the same for all times. 

Next we study the constant Cd- First we use critical 
temperature initial conditions using the dynamic results 
to avoid the ambiguity introduced by the units restoring 
constant in equilibrium. In short we find 



Note that this value is slightly higher than Ch, in accord 
with the prediction from the sum rules and the analytic 
argument, and the difference between the two {cd — Ch ^ 
0.002) is slightly smaller than the one that is obtained 
from the sum rules. 



C. Effect of the working temperature 

Up to now we have considered quenches to a zero work- 
ing temperature. In this Section we investigate the effect 
of having a finite temperature on the dynamics. 

The arguments in Sect. IIVI rely on the T = AUen- 
Cahn equation (|29p . Temperature fiuctuations have a 
two- fold effect. On the one hand they generate equilib- 
rium thermal domains that arc not related to the coars- 
ening process. On the other hand they roughen the do- 
main walls thus opposing the curvature driven growth 
and slowing it down. 

Renormalization group treatments of domain growth 
dynamics (23| have led to the idea that a T = fix point 
controls the domain growth for all T < Tc, i.e. that 
thermal fiuctuations are irrelevant to the asymptotic dy- 
namics of the ordering system, their contribution being 
limited primarily to the renormalization of temperature- 
dependent prefactors. 

For the distribution of domain areas and hull enclosed 
areas, one may expect that once equilibrium thermal do- 
mains are subtracted - hulls and domains associated to 
the coarsening process are correctly identified - the full 
temperature dependence enters only through the values 
of Xh and A^, which set the time scale. 

The first step then is to identify the temperature de- 
pendence of the parameter A^- The simplest and most 
direct way to do this is to use the scaling hypothesis and 
analyze the behaviour of the spatial correlation 

1 ^ 

C{r,t) = ^^(Si(Osj(i))||n-r,|=r 



R{t)J ' 



(59) 



Cd 



0.0251 ± 0.0015 



(58) 



where m(T) is the equilibrium magnetization density and 
a < r < L and to < i- Using R{t) - [A<j(T)t]i/2, the 
T-dependence of Xd can be estimated either by collaps- 
ing all curves or by studying the value of r at which 
C{r,t) = 1/2. The resulting Xd{T) obtained using these 
two prescriptions is shown in Fig. 1171 Arf(T) is a mono- 
tonically decreasing function of temperature, starting at 
Xd{T = 0) = 2.1 and falling-off to zero at T^. These re- 
sults are consistent with the evaluation of Xd.h{T) from 
the analysis of nd^h{A,t), see below, and it is at variance 
with what was previously presented in Refs. [2^ . 

Assuming that Xd vanishes at Tc one can deriye the 
way in which it does with a simple argument [2^. We 
require that the coarsening law for coarsening below 
namely R{t) - [X{T)t]^^^, match critical coarsening at 
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FIG. 17: (Colour online.) Top panel: spatial decay of the 
equal-time correlation, Eq. (|59p . at fixed T — and several 
diff'erent times t. Bottom panel: the T dependence of the 
parameters Xd,h- Two sets of data points are extracted from 
the analysis of the correlations shown in the top panel. The 
data named half- value are obtained from C{r,t) = 1/2 and 
the data named collapse from collapsing the curves on the 
range r ^ a. The two other sets are obtained from fitting 
nd,hiA,t) with Ad,/i(T) as free parameters. 

Tc, viz. R{t) ^ i^/^ with z the dynamic exponent, for 
T Tc- Near (but just below) Tc the coarsening length 
grows as ^-''{T)t^^^ as long as R{t) > £_{T) with ^(T) 
the equilibrium correlation length. For R{t) comparable 
with £,{T), this has to be modified by a function of R{t)/£^ 
and, since R{t) ~ t^^^ at Tc, we can write 

i?(0~r"(r)iV2/(^_l_^ . (60) 

In the limit ^{T) oo, the £, - dependence must drop out. 
In order to cancel the time dependence at large times, one 
needs f{x) ~ x^/^ for a; ^ oo. This yields R{t) ~ f/^, 
which fixes the exponent 'a' as a = (2 — z)/2, giving 
A(r) - ('^''{T) - {Tc - Ty^^''^\ Inserting the exact 
value v = 1 and the numerical value z — 2.15(2) p6l | 
implies 

Xd{T) ^ {Tc - r)"-i5 . (61) 



FIG. 18: (Colour online.) The number density of hull enclosed 
areas after t — 128 MCs at the working temperatures T = 
0.5, 1.5, 2 and T = T^. 

Note that we are matching two nonequilibrium growth 
laws - the one below Tc and the one at Tc - not an equi- 
librium and a nonequlibrium one. The data in Fig. 1 171 are 
still far from the critical region where this small power- 
law decay should show up. 

Finite working temperatures also affect the distribu- 
tion of domain areas. In Fig. [TH] the raw data at i = 128 
MCs is shown for four working temperatures. Upward 
deviations with respect to the result of zero working tem- 
perature are prominent in the small areas region of the 
figure, and increase with temperature. 

In Fig. [Tni (top) we display the raw data at the work- 
ing temperature T = 1.5, for several times. Notice that 
although the curves move downwards, the small areas re- 
gion becomes time independent. This region also fails 
to collapse (bottom) with the proposed scaling using 
the temperature dependent values of Xd{T). The reason 
is that the distribution counts thermal equilibrium do- 
mains, that is to say fluctuations that are present in an 
equilibrated sample at the working temperature, but are 
not due to the coarsening process. Thus, these fluctua- 
tions should be identified and eliminated from the statis- 
tics. We tried to apply the method introduced by Der- 
rida [5| , and extended by Hinrichsen and Antoni [23| , to 
eliminate thermal domains, but the results were not sat- 
isfactory, as not all of them could be eliminated. Thus, 
instead of removing each thermal domain, we tried to di- 
rectly remove their contribution to the distributions by 
simulating samples in equilibrium at the working tem- 
perature, starting with a fully magnetized state, and 
computing the number density of thermal domain ar- 
eas. These data are shown with green data points in 
Fig. [19] and [20l Surprisingly enough, thermal fluctua- 
tions generate areas that are larger than one would have 
naively expected. Equilibrium arguments suggest that 
the averaged area of thermally generated domains scale 
as At_~ e{T) with ^{T)/po ~ /-(I - T/Tc)-", v - 1 
and = 0.18 [28]. This estimate yields, for example. 
At ^ 4ylo at T = 1.5. In equilibrium at this tempera- 
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FIG. 19: (Colour online.) The number density of hulls for 
T = 1.5 after different times (top) and the scaling of these 
data points (bottom) 



ture the average size of the domains found numerically is 
(At) ^ 1.5Ao. However, the probability distribution of 
thermal areas has a non-negligible weight - as compared 
to the one of coarsening domains ~ that goes well beyond 
this value. For example, in Figs. fT9l and [20l we see that 
the crossover between the thermal area distribution and 
the coarsening area distribution occurs at A ~ 10 Aq. 

In Fig. [201 we also present data for the dynamic distri- 
bution at three different times, and compared with the 
analytic prediction using Aci(T) estimated from the anal- 
ysis of the global spatial correlation, see Fig. [T71 We 
conclude that the agreement between analytic prediction 
and numerical results is very good in the region in which 
the thermal domains are subdominant, i.e. when the 
blue dynamic curves deviate from the green equilibrium 
one, indeed the regime in which the analytic calculation 
is expected to apply. 

One can also use the results in Fig. [21] to estimate 
the value of Ad(r). Indeed, a fit of the numerical data 
for areas larger than the value at which the equilibrium 
thermal contribution (green points) deviates from the dy- 
namic one, yields the values of Xd{T) [and A/i(r)] shown 
in Fig.[T7l This analysis allows us to extract independent 
predictions for Xd{T) and Xfi{T). We find that the quali- 
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FIG. 20: (Colour online.) The contribution of 'thermal do- 
mains' obtained by simulating an equilibrated sample at the 
working temperature T — 1.5, along with the evolution of 
the distribution of domain sizes after a quench to the same 
temperature. 
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FIG. 21: (Colour online.) Zoom on the number density of 
domain areas at t = 128 MCs at three working temperatures 
given in the key. The black lines are the equilibrium distri- 
butions at r = 1.5 and T = 2 and the other lines (pink, blue 
and green) represent our analytic prediction for the coarsening 
areas. 



tative T-dependence is the same. As regards the absolute 
values, the numerical data yield Xh{T) < Xd{T) on the 
whole range. Note that the sum rules suggested A^; — Xh 
and the analytic prediction Xd = Xh + 0{ch). 



VI. STATISTICS OF PERIMETERS AND 
FRACTAL PROPERTIES 

The analytic argument described in Sect. [TV] can be ex- 
tended to study the distribution of domain wall lengths 
or perimeters. In this Section we present the analytic 
prediction for this function together with numeric results 
that confirm it. We study two types of domain bound- 
aries: those associated to the hulls and those associated 
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to the domains that is to say that include external and in- 
ternal perimeters. In the simulations we define the length 
of the boundary as the number of broken bonds. 



A. Initial conditions 

1. Equilibrium at Tq = Tc 

In equilibrium wc find numerically that the domain 
areas and their corresponding boundaries are related by 
(see Fig. [22] where the scatter plots have been averaged 
to make the trend clearer) 

Ah cl'^ with al'^ - 1.47 ±0.1, 



A, 



in the whole range of variation. Note that the longest 
lengths, p w 10^ — 10** may be affected by finite size 
effects given that the linear size of the simulating box is 
L = 10^. The spanning clusters are not counted (note 
that their perimeters would be severely under estimated 
due to the periodic boundary conditions). The exponent 
•^d ~ 1.47 ± 0.1 is consistent with the result in [l3| 
mentioned in Sect. Ill CI The constants take the values 
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between the small and large p regimes in the relation 
between areas and perimeters for the domains is due to 
the existance of holes in the large structures. The small 
domains and hulls are just the same objects because the 
former do not have holes within. 
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FIG. 22: (Colour online.) Relation between areas and perime- 
ters in equilibrium at To = Tc- 



Numerically, we find that the number densities of hull 
and domain lengths at critical Ising conditions are (see 
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FIG. 23; (Colour online.) Distribution of domain and hull 
lengths in equilibrium at To — Tc. 



Fig. 

nhip,0) - p-^'^ with a --2.48 ±0.05. 

p-C , (> ^ 2.17 ± 0.05 for p ~ 50 , 



nd{p, 0) 



p-C , (< ^ 2.48 ± 0.05 for p - 50 



The value of is to be compared to the analytic result 
= 27/11 « 2.454 [l3|. It is interesting to notice that 
the distribution of domain lengths is not a single power 
law in constrast to the distribution of domain areas. 



S. Equilibrium at To — oo 

After a few time-steps evolving at T = from the 
infinite temperature initial condition, we reach critical 
percolation conditions. 

In Fig. I24| we show the area perimeter relation for 
hulls, in equilibrium at To — > oo and after a few time 
steps. The analysis of this figure and the corresponding 
one for domains yields 



Ah 



Ad 



with 



P " , 



1.12 ±0.1 



1.01 ±0.1 for p- 50 , 



for p ^ 50 , 



The constants take the values cl = 0.96, and c^* '^^ 



h — O'ii'-t '^d 

1.50. The exponent a); cannot be determined numer- 
ically since critical percolation is not accessible exactly. 

As we can see in Fig. [55l the initial weight of the num- 
ber density at large values of the perimeter is lower than 
expected at critical percolation. In a few time-steps long 
perimeters develop and the weight reaches the asymp- 
totic power law at large values of p while it looses weight 
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FIG. 24: (Colour online.) Relation between areas and perime- 
ters in equilibrium at To — > oo and after a few time steps when 
critical percolation is approximatively reached. 




P 

FIG. 25: (Colour online.) The number density of hull lengths 
in equilibrium at To — > oo labelled as t = and the same 
quantity evaluated at two later times after evolution at zero 
working temperature. The lines are nh ~ p~^'^^\ see the text 
for a discussion. 



at small values of p. This effect is the same as the one 
observed in the study of the initial and early times num- 
ber density of areas, see the discussion in Sect. IV Al and 
Fig. [5] For the hull and domain length distributions one 
finds 

nft(p,0) - p-^'>^ with 2.12 ±0.05 

{p-^'-i^ , C> - 2.01 ± 0.05 for p - 50 , 
p-'^i , Cd for p - 50 , 

The analytical result for the hull exponent in critical per- 
colation is = 15/7 « 2.14 [TT^. 



3. General comments on both initial cases 

It is interesting to note that the exponents character- 
ising the number density of perimeter lengths at the two 
initial conditions are significantly different. They are ap- 
proximately equal to 2.5 at Tc and 2 at Tq — s- oo. This is 
to be contrasted with the behaviour of the area number 
densities for which the exponents were identical for hull 
enclosed areas and very close indeed for domains. 

The exponents a and C are linked by the fact that each 
hull-enclosed area or domain area is in one-to-one relation 
to its own boundary. Thus, nh{A, Q)dA = nh{p, 0)dp and 
one finds 

n,(p,0)^^p-i-"i", (62) 

which implies 

a = 1 + c/^^ , (63) 

These conditions are also satisfied for the primed (To 
oo) quantities. Within our numerical accuracy these re- 
lations are respected, for instance 

~ 2.12 , ^ 1.12 , To ^ cx) , 

a ~ 2.48 , ai*^ - 1.47 , To = T, 

Similarly, for domain areas and domain boundaries one 
obtains 

n4P,0) - ^ p-'-^^-'^-'^\ (64) 

therefore, 

Crf - 1 + (r - l)a« , (65) 

These relations are satisfied for both and as 

well as for primed [Tq — » oo) quantities. They are re- 
spected by our measures. 

The main sources of error in the determination of the 
exponents and the constants in the study of the initial 
conditions are the following: (i) statistical errors, al- 
though we have a rather good sampling; (ii) the choice 
of the large area-perimeter limit that is not perturbed 
by finite size effects, and (iii) the fact that the Tq ^ oo 
initial condition is not exactly at critical percolation. We 
estimated the magnitude of the error to be ±0.1 in the 
a exponents, and ±0.05 in the C exponents, which corre- 
spond to, roughly, less than 10% in both cases. Within 
this level of accuracy, the relations between exponents 
([S5)) and are satisfied. 

B. Time evolution at zero temperature 
1. Hulls 

After a quench from Tq = Tc, the hull enclosed areas 
and their corresponding perimeters, during coarsening at 
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zero temperature, obey the scaling relations (see Fig. 
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FIG. 26: (Colour online.) Time evolution of the hull enclosed 
area vs. perimeter relation for To = Tc and different times 
indicated in the legend. 




with 



and 



- 1.37±0.2 
7j> ~ 0.26 



a< - 1.83 ±0.2 
r/< - 0.06 



for - 50 , 

Xht 



for A < 10 , 
Xht 



and, after a quench from Tq — s- cx), 



with 




a'>^ 1.12 ±0.2 1 ^ > 

V't - 0.38 / ^"'^ A;:^ 



(66) 
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FIG. 27: (Colour online). Hull enclosed area vs. perimeter 
relation at time £ = 32 MCs from two different initial condi- 
tions. 



and 



1.83 ±0.2 



rj'< ^ 0.057 



for 



A 
X^ 



10 



(71) 



We note that the relation between area and perimeter 
exhibits two distinct regimes. During the coarsening pro- 
cess a characteristic scale A* {t) ~ Xht develops such that 
domains with area A > A* have the same exponent as in 
the initial condition (structures that are highly ramified 
with a smaller than two) and domains with A < A* are 
regular (a ~ 2) (as shown in Fig.[27]the structure of these 
small domains does not depend on the initial condition). 
This phenomenon is reminiscent of an unroughening tran- 
sition occurring at a velocity . The same features were 
observed by Grest and Srolovitz [2^ and Fialkowski and 
Holyst [30] in the study of the domain fractal dimension 
during coarsening. 

The hull structures of any size do not have holes, there- 
fore the crossover we see is of pure dynamical origin, con- 
trasting the idea presented in j3l| , where the crossover in 
the domains where explained by only geometrical rea- 
sons. 

Note that we estimated the error in the exponents a 
to be ±0.2 and thus more important than in the analysis 
of the initial conditions. The reason is that the crossover 
from the small area to the large area regime is not suf- 
ficiently sharp and the choice of the fitting interval in- 
troduces an additional source of error. Indeed, note that 
in Eqs. (pS)) and we did not use the intermediate 
regime 10 < ^ < 50 to fit the power laws. 

In analogy with the derivation in Sect. lIVl for the time- 
dependent number density of domain areas, the time- 
dependent number densities of hull and domain wall 
lengths are given by 



nh,d{p,t) ^ dp^ 6{p- p{t,pi))nh,d{Pi,ti) (72) 
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with nh,d(pi,ti) the initial condition and p{t,pi) the 
perimeter length of a boundary at time t that had initial 
length Pi at time U . 

Let us here discuss the hull lengths. In this case one 
can simply use the exact number density of hull enclosed 
areas, nh{A,t) ^ Ch/{A + Xht)^ for, say, To = Tc and 
Eq. (j66p to relate time-dependent areas to their perime- 
ters on the two regimes of large and small areas. After a 
little algebra one derives 



(Xhtf/^ Uhip, t) 



(73) 



for small areas, A/Xht < 10, and 



( V 



l + r/> 



(74) 



for large areas AjXht > 50. Note that these expressions 
satisfy scaling - see Eq. (fTTj) . Interestingly, the scaling 
function, /<(a;) = a;"^~^/(l-|-?7^a;"^ )^ with x = p/^/Xhi 
reaches a maximum at 



1 



af + 1] 



(75) 



and then falls-off to zero as another power-law. There 
is then a maximum at a finite and positive value of p 
as long as > 1, that is to say, in the regime of not 
too large areas. The numeric evaluation of the right- 
hand-side yields Xmax = Pmax / [s/ Xht) ^ 3 which is in 
the range of validity of the scaling function /< . The be- 
haviour of the time-dependent perimeter number density 
for long perimeters is controlled by Eq. ([71)1 that falls-off 
as a power law />(x) ^ x^^^^"-^) , Although the function 
/> also has a maximum, this one falls out of its range of 
validity. 

Above we used the critical Ising parameters. The re- 
sults after a quench from Tq ^ oo follow the same func- 
tional form with the corresponding primed values of a 
and 77 and Ch '2cfi. 

The power law describing the tail of the number den- 
sity of long perimeters is the same as the one character- 
ising the initial distribution, since = aj^*' and then 
1 + a-f^ — Ch- Therefore, the decay of the time-dependent 
number density at long perimeters after a quench from 
To = Tc and Tq 00 are distinguishably different with 

~ 2.5 and « 2. This is to be contrasted with the 
small difference in the area number densities that fall 
with two power laws that are so close (powers of 2 and 
2.05) that are impossible to distinguish numerically. 

In Fig. [251 top and bottom, we display the time- 
dependent perimeter number densities for a system evolv- 
ing at zero temperature after a quench from Tq = Tc and 



To — > 00, respectively. Notice that the perimeter length 
definition we are using on the lattice can only take even 
values and thus when constructing the histogram we have 
to take into account the extra factor of 2 in the binning. 
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FIG. 28: (Colour online.) The time-dependent number den- 
sity of perimeters evolving at T = from an initial condition 
at To — Tc (top) and To — > 00 (bottom). Note that the time- 
dependence is visible in the whole range of values of p (while 
in the area number densities the large area tails were very 
weakly dependent on time, see Fig. [7}. 

In Fig. [211 we display the scaling plot of the number 
density of hull lengths and we compare it to the analytic 
prediction ([73| and ([74)1 . The data are in remarkably 
good agreement with the analytic prediction; the lines 
represent the theoretical functional forms for long and 
short lengths, and describe very well the two limiting 
wings of the number density. The maximum is located at 
a value that is in agreement with the prediction, Eq. ([75)1 . 



2. Domains 

We studied the relation between domain areas and 
their corresponding perimeters during coarsening at zero 
temperature finding that the scaling forms 



A 

x7t 
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(76) 



22 





FIG. 29: (Colour online.) Scaling of the time-dependent num- 
ber density of hull lengths evolving at T = from an initial 
condition at To = Tc (top) and Tq ^ oo (bottom). The 
solid black lines represent the theoretical prediction valid for 
A/\ht < 10 and for A/\ht > 50. The agreement between the- 
ory and numerical data is again very impressive. The small 
grey line in the top plot represents the slope in the bottom 
plot, showing that in contrast with domain size distribution, 
perimeter distribution is very sensitive to initial conditions. 
The isolated data points that lie above the scaling function 
correspond to reversed, isolated spins within a bulk of the op- 
posite sign that give rise to a perimeter length p — 4 (four 
broken bonds). The area number densities also showed this 
anomalous behaviour for A = 1. 



with 




FIG. 30: (Colour online.) The time-dependent relation be- 
tween the area and the perimeter of a cluster of aligned spins 
(domain) evolving at T = after a quench from To = Tc (top) 
and To — > oo (bottom). In both cases Ad = 2.1. 



and 



^^^^ These results are shown in Fig. [301 Note that the large 

Q,> I IQ ^ j[ ^^^^ results match the behaviour of the initial conditions 

> n r ^ ' (''"'') both cases, and small domains are much more compact 

With the same line of argument exposed above we can 
analyse the statistics of the domain walls, that is to say, 
^ 1 83 1 A including external and internal perimeters. One finds 

^ n nc;7 f ^'^^ T~r ^ ' C''^) basically the same results as for the hulls; for critical 

Ising initial conditions: 



a 



^ 0.057 ; X,t 
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for small areas and its obvious modification for large ar- 
eas. For To ^ oo one replaces 77^ and ad by the primed 
quantities and Cd — > 2cd. The scaling analysis of the 
number density of domain wall lengths is displayed in 
Fig. [31] for both initial conditions. Once again we find 
a very good agreement between the analytic predictions 
and the numerical data. 



exponent (Tq ^ 00 or Tq = Tc). This is reasonable since 
the large structures are still 'unaware' of the coarsening 
process and thus retain the form they had in the initial 
configuration. For small areas, instead, we see domain 
walls roughening due to thermal agitation but it is hard 
to extract the value of the exponent with sufficient 
accuracy. 



3 



->3 



"a 




FIG. 31: (Colour online.) (Colour online.) Scaling of the 
time-dependent number density of domain wall lengths evolv- 
ing at T = from an initial condition at To = Tc (top) and 
To — » 00 (bottom). The solid black lines represent the theo- 
retical prediction valid for A/\ht < 10 and for A/X^t > 50. 
The origin of the isolated data points is the same as in Fig. [29] 



C. Finite temperature evolution 

Once we analysed the statistics of perimeters in the 
zero temperature dynamics we focus on the effects of a 
finite working temperature. We briefly list the results 
below without presenting the data. 



1. The area-perimeter relations 

For large areas we find the same exponent as for zero 
temperature coarsening that is also the initial condition 



2. Perimeter number densities 

The scaling of the perimeter number densities and the 
functional form for the scaling function predicted ana- 
lytically describe the numerical data with high precision 
once the values of the exponents a, the prefactors rj, 
and the parameter A are modified to take into account 
thermal agitation (for absolute area values larger than 
A ~ lO^Oi this value limits the range of areas where the 
effect of thermal fluctuations is larger than the ones of 
coarsening). The analytic prediction is very accurate in 
the region of small coarsening domains, A/Xd^ht < 10 
and A > IOAq where the maximum is located, and in 
the region of large coarsening domains, A/Xd,ht > 10 
and A > lO^o, for both domains and hulls and the two 
initial conditions. 



VII. CONCLUSIONS 

In this paper we studied the statistics and geometry of 
hull enclosed and domain areas and interfaces during the 
non-equilibrium dynamics of curvature driven pure coars- 
ening in two dimensions. The analytical part of our work 
relies on the AUen-Cahn equation derived from the con- 
tinuous Ginzburg-Landau field-theory in two-dimensions 
while the numerical part of it dealts with Monte Carlo 
simulations of the 2c?IM. Our main results are: 

(i) We proved scaling of the various number densities 
studied. 

(ii) We derived the exact number density of hull enclosed 
areas and hull lengths; we obtained approximate expres- 
sions for the number density of domain areas and domain 
wall lengths. 

(iii) The geometrical properties and distribution of the 
time dependent large structures (by large we mean much 
larger than the average ones) are the ones of critical con- 
tinuous percolation (for all initial conditions equilibrated 
at To > Tc) and critical Ising (for Tq = Tc). The long in- 
terfaces retain the fractal geometry imposed by the equi- 
librium initial condition and the scaling function of all 
number densities decay as power laws. 

(iv) Instead, small structures progressively become regu- 
lar and the area-perimeter relation is A ~ . 

(v) We took into account the effects of a finite work- 
ing temperature by correctly eliminating purely thermal 
fluctuations and thus correctly identifying the coarsening 
structures. The temperature effect thus amounts to in- 
troducing the temperature dependence in the prefactor 
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in the growth law, R{t) - [X{T)t]^/'^. A(T) is a monoton- 
ically decreasing function of T that vanishes at T^. 

It is important to stress that our analytic results rely 
on the use of the AUen-Cahn result for the velocity of an 
almost flat interface. Thus, they would be expected to 
hold only in a statistical sense and for large structures in 
the lattice model. Surprisingly, we found with numerical 
simulations that the number density area distributions in 
the 2dlM match the analytic predictions for very small 
structures, and even after a few MC steps evolution of 
a critical Ising initial condition for which rather rough 
interfaces exist. 

Using the AUen-Cahn result and a variety of numerical 
measurements we verified the well-known result for pure 
coarsening with non-conserved order parameter: there is 
a characteristic growing length that increases in time as 
i^/^. The mesoscopic analysis presented here allows us to 
demonstrate that the reason for the growth of the char- 
acteristic length is the disappearance of small structures. 

Our analytic results hold only in two dimensions. In 
Appendix [B] we summarized the behaviour of the dis- 
tribution of domain lengths in the Ising chain. As ex- 
pected from scaling arguments, this quantity scales with 
the typical domain length, R{t) ~ t^^^ but the form of 
the scaling function is very different from the one in two 
dimensions: the pdf vanishes at zero scaling argument 
{x = 0), it then increases linearly to reach a maximum 
and then falls-off to zero exponentially. On the other 
hand, we cannot extend the analytic argument to dimen- 
sions higher than two since the hull enclosed volumes no 
longer decrease with time in a manner that does not de- 
pend on their own volume (20| . 

This paper and [J] open the way to a number of re- 
lated studies. For instance, it would be interesting to 
extend the analysis presented here to dynamic clusters 
of correlated spins (droplets), that are known to describe 
thermal fluctuations close to the transition (see [3§| and 
pol l for a review). These droplets are smaller than the 
geometric clusters in which they are embedded because 
some of the neighboring parallel spins are discarded (by 
a temperature dependent criteria) for not being corre- 
lated (remember that even at infinite temperature, where 
all spins are uncorrelated, there are domains of parallel 
spins). 

Two-dimensional coarsening with conserved order pa- 
rameter is another problem that deserves a careful study 
along these lines. 

In p[ we include a preliminary analysis of the hull 
and domain structure in the finite temperature dynamics 
of the bi-dimensional random bond ferromagnetic Ising 
model after a quench from infinite temperature. In such 
a disordered case a finite working temperature is nec- 
essary to help the interfaces depin from pinning cen- 
ters in the quenched disordered potential through ther- 
mal activation. We first computed the typical domain 
radius that scales the time-dependent spatial correla- 
tion, C{r,t) ~ f{r/R{t)), in the scaling regime. Due 
to the presence of quenched disorder R{t) strongly de- 



pends on temperature and the strength of randomness 
and it is slower than the simple square root behaviour 
of the pure Ising case. We then showed numerically that 
the number density of hull enclosed and domain areas 
scale as R'^{t)nh.d{A,t) ~ g{A/ R'^{t)) for areas satis- 
fying 10-1 - A/R^{t) and that are smaller than the 
cut-off set by finite size effects. The effect of a non- 
trivial typical radius R{t) determined by the quenched 
disorder can be tested in the intermediate regime, say 
10^1 ^ A/R^{t) ^ 10^, where the quality of the scahng 
plot is excellent. The scaling function g{x) does not de- 
pend on the disorder strength satisfying the hyper-scaling 

hypothesis [if. For smaller areas, say A/R^{t) ^ 10-1, 
the contribution of thermal domains with domain walls 
roughened by disorder is important. We shall give more 
details on the domain morphology of the quenched disor- 
dered coarsening problem in a separate publication 42] . 

These results give an idea of the richness and com- 
plexity of coarsening phenomena even in the absence of 
quenched randomness. We expect them to be of help 
in understanding the fluctuating dynamics of even more 
complex situations, like spin-glasses and glassy prob- 
lems (4H , in which the mere existence of a domain growth 
of two competing equilibrium phases is not even estab- 
lished. 



APPENDIX A: ALGORITHM USED TO 
IDENTIFY HULL ENCLOSED AREAS 

In order to obtain the size of each hull in an L x L 
system, a biased walk along the interior border of each 
domain is performed, with the hull enclosed area being 
updated at each step. This algorithm is related to, but 
different from the Grossman- Aharony proposal. 

Labelling. The N sites are initially indexed from 
to A'' — 1 (top-left=0) while the domains (geometric, 
Coniglio-Klein, etc) are identified and labelled by the 
Hoshen-Kopelman algorithm [21|]. By construction, all 
sites in each domain receive the (unique) label corre- 
sponding to the smaller index among its spins. 

Starting point. The putative starting site for the walk 
is the spin whose index identifies the cluster. In some 
cases (for example, when the cluster crosses a border), 
it may not be the leftmost /topmost site, as it should 
in order to be counted correctly by the implementation 
of the algorithm below. Although such domains may be 
excluded from the statistics, for finite sizes the introduced 
bias is unacceptable. Thus, before starting, we try to find 
another site in the same cluster above or to the left of the 
original spin, where a new starting site may be found. 
Once this first site is correctly identified, we assign a 
height yo- 

The walk. From the starting point, we try to turn 
clockwise around the domain border. Viewed from the 
incoming direction, the attempted move is performed in 
the sequence: left, front, right and backwards. For the 
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TABLE I: (Colour online.) Incremental area contribution for 
each step during the oriented walk along the domain border. 
The values depend both on the present (t) and former steps 
(t—l). On the right we show the labelling of the four possible 
directions. 

sake of notation, we label the four directions with the 
indices shown in the table [XI Thus, for example, if the 
previous step was 0, the attempted order will be 1, 0, 3 
and 2. 

The area. The first step sets the area to the value 
^1 = yo + 1- As the walk proceeds, both and area and the 
height are updated, Ai^i = Ai + AA and yi+i = yi + Ay 
(in this order), where and Ay depend also on the 
former direction (see table E]) . At the end of the walk, 
when the departure site in reached, if the last direction is 
1, we increase A by —y. Care should be taken when the 
starting point has right and bottom neighbours belonging 
to the same cluster (but not the bottom right diagonal), 
because in this case the walk should only be finished after 
the second time the starting point is visited. 

APPENDIX B: THE ONE DIMENSIONAL CASE 

For the sake of comparison we present in this Ap- 
pendix results for the domain size distribution in the 
one-dimensional Ising chain. The Id ferromagnet with 
nearest neighbours interactions is given by the Hamilto- 
nian 

W = -J^s,s,+i (Bl) 

i 

where the spins are Si — ±1 and J > 0. The system 
evolves through a Heat Bath dynamics, where each spin 
tries to align with its local field (with a temperature de- 
pendent probability). 
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FIG. 32: (Colour online.) The average domain length at finite 
temperature, = lyfwi independently of temperature in 
the domain- growth regime before saturation. 
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FIG. 33: (Colour online.) Collapsed distributions of the 
(rescaled) probability distribution of domain lengths after a 
quench to T = for L = 10^ and times t = 2^, . . . , 2^. AU 
data collapse onto an universal curve whose tail is exponen- 
tial 36,]. Inset: the distribution for short length scales where 
the universal function is j(x) ~ -kx. Notice that the average 
domain length is different from the typical one. 

We measure the domain length distribution at time 
n(^, t) (\ < I < V) after quenching the system from 
To — > oo to a final, working temperature, T. There is 
no finite static critical temperature and the system or- 
ders ferromagnetically only at T = 0. At all finite tem- 
peratures a completely disordered initial condition starts 
evolving in a coarsening regime in which regions of fi- 
nite length order and subsequently crosses over to equi- 
librium in the paramagnetic phase. Differently from the 
higher dimensional case where domain growth is driven 
by interfacial tension, in c? = 1 coarsening is driven by 
the diffusion of domain walls and annihilation when they 
meet. 

In the infinite temperature initial condition the spins 
are uncorrelated and n(^, 0) — 2^^, while in equilib- 
rium, the normalized distribution is also exponential 
as it corresponds to the distribution of domains of a 
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FIG. 34: (Colour online.) Domain length distribution at T = 
0.5 for increasing values of time (from left to right; 0, 4, 16, 
64, 256 and oo. The initial, n{£,0), and final distributions, 
n{£,oo), are also indicated with solid lines (see text). The 
small £ region is shown in the inset. 



paramagnet in a field: n{i, oo) = r(l — r) , with 
r= [1 + cxp(2/3)]-i p. 

During the coarsening regime, the distribution of do- 
main sizes obeys the scaling behavior 



im) 



(B2) 



where the time dependent average length is {£{t)) ~ 
2\pKi [s^, [SI], independently of temperature and satu- 
rates at {£{oo)) = 1/r after a time that roughly grows as 
exp(4/T), see Fig. [321 At zero working ternperature the 
universal function f{x) is given by 34 [36L IstI. Issjl 




-Ax + B) 



a; < 1 , 
a; > 1 , 
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where the constants are known exactly: A = 1.3062 and 
B = 0.597 [ail, see Fig. [3S Differently from higher di- 
mensions, here as large clusters coalesce, there is no for- 
mation of small clusters (compared with the typical clus- 
ter). This fact is at the origin of the dip in the pdf close 
to = 0. 

At finite working temperatures, because of the thermal 
fluctuations, small clusters are created at a constant tem- 
perature dependent rate. As say one spin flips within a 
domain not only there is a new domain with length 1=1 
but the host domain has been cut in two pieces of relative 



much shorter length. Thermal agitation thus decreases 
the depth of the dip. When x > 1, f{x) is still well fit- 
ted by Eq. (jB3p . In the opposite regime, x ^ 1, instead 
thermal fluctuations diminish the dip and the behavior 
deviates from the linear scaling function, as can be seen 
in Fig. El 
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FIG. 35: (Colour online.) Temperature dependence of the 
domain length distribution at small lengths. Filled symbols 
are for T = Q while hollow ones are for T = 0.5. The line is 
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